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Abstract 

The  presented  exact  solutions  of  homogeneous  Maxwell’s  equations  in  complex  media  and 
free  space  describe  fields  having  a rather  involved  curl  structure  and  a very  small — about 
several  wave  lengthes  of  composing  plane  waves — and  clearly  defined  core  region  with  max- 
imum intensity  of  field  oscillations.  In  a given  Lorentz  frame  L,  a set  of  the  obtained  exact 
time-harmonic  solutions  of  the  free-space  homogeneous  Maxwell  equations  consists  of  three 
subsets — termed  “storms” , “whirls” , and  “tornadoes”  for  the  sake  of  brevity — for  which  time 
average  energy  flux  is  identically  zero  at  all  points,  azimuthal,  and  spiral,  respectively.  In 
any  other  Lorentz  frame  L',  they  will  be  observed  as  a kind  of  electromagnetic  missiles  mov- 
ing without  dispersing  at  speed  V < c.  The  solutions  which  describe  finite-energy  evolving 
electromagnetic  storms,  whirls,  tornadoes  are  also  presented. 


1.  Introduction 

In  the  beginning  of  eighties,  Brittingham  [1]  proposed  the  problem  of  searching  for  specific 
electromagnetic  waves — focus  wave  modes — having  a three-dimensional  pulse  structure,  being 
nondispersive  for  all  time,  and  moving  at  light  velocity  in  straight  lines.  A number  of  packet- 
like solutions  have  been  presented  [1,  2,  3].  However,  it  seems  [3,  4,  5],  finite-energy  focus  wave 
modes  can  not  exist  without  sources.  In  1985,  Wu  introduced  [5]  a conception  of  electromagnetic 
missiles  moving  at  light  velocity  and  having  a very  slow  rate  of  decrease  with  distance.  Because 
of  these  properties,  such  missiles  have  important  possible  applications  [5]. 

In  Ref.  [6],  vector  plane-wave  superpositions  defined  by  a given  set  of  orthonormal  scalar 
functions  on  a two-dimensional  or  three-dimensional  manifold — beam  manifold  B — are  treated. 
The  proposed  approach  makes  it  possible  to  compose  a set  of  orthonormal  beams,  normalized  to 
either  the  energy  flux  through  a given  plane  <jq  with  unit  normal  q (beams  with  two-dimensional 
B)  or  to  the  total  energy  transmitted  through  this  plane  (beams  with  three-dimensional  B),  as 
well  as  some  other  specific  exact  solutions  of  wave  equations  such  as  three-dimensional  standing 
waves,  moving  and  evolving  whirls.  This  approach  can  be  applied  to  any  linear  field,  such  as 
electromagnetic  waves  in  free  space,  isotropic,  anisotropic,  and  bianisotropic  media,  elastic  waves 
in  isotropic  and  anisotropic  media,  sound  waves,  etc.  By  way  of  illustration,  some  specific  families 
of  exact  solutions  of  the  homogeneous  Maxwell  equations,  describing  localized  electromagnetic 
fields  in  free  space,  are  obtained  in  Refs.  [6,  7,  8].  In  this  paper,  we  present  some  new  types  of 
such  localized  electromagnetic  fields  in  complex  media  and  free  space. 
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2.  Time-Harmonic  Localized  Fields  Defined  by  the  Spherical  Harmonics 

In  this  paper,  we  treat  time-harmonic  electromagnetic  fields  in  linear  complex  media  or  free 
space,  defined  by  the  spherical  harmonics  as 

/■27T  /*$2 

W2(r,t)  =e~iut  / dtp  eiT-k^Y?(d,tp)v(e,tp)W(0,<p)sm0d9,  (1) 

Jo  Je  i 

where  W = col(E,  B)  specifies  the  polarization  of  the  eigenwave  with  the  wave  vector  k.  To 
compose  these  beams,  it  is  necessary  first  to  calculate  parameters  of  eigenwaves.  The  corre- 
sponding relations  for  electromagnetic  waves  in  a general  bianisotropic  medium  are  presented  in 
Ref.  [6]. 

There  are  two  main  ways  to  set  the  beam  base,  i.e.,  to  specify  the  functions  k = k(0,  tp) 
and  W = W(0,  tp).  One  can  set  first  the  unit  wave  normals  of  these  eigenwaves  by  a function 
k = k(0,  tp).  In  particular,  one  can  set  the  angular  spectrum  of  plane  waves  by 

k = k/fc  = sin0'(ei  cos  tp  + e?  sin<^)  + e3  cos  9 , (2) 

where  9'  = kq9,  and  kq  is  some  real  parameter.  Then,  one  has  to  calculate  the  refractive  indices 
nj(9 , p)  = nj(k(9,  p))  of  all  isonormal  waves  and,  by  choosing  some  branch  rij(9,  p),  to  specify 
the  wave  vector  function  k(6)  = (cu/c)rij(9,  p)k(9,  p)  and  the  amplitude  function  W(0,p)  = 
col(E(0,  p),  B(0,  p))  as  well.  The  alternative  is  to  set  first  the  tangenial  components  of  wave 
vectors  by  a real  vector  function  t = t(0,  p).  Then,  the  normal  component  £j(0,  p)  — £j( t(0,  p)) 
of  k(0,  p)  — t (9,  p)  + £j{9,  p)*i  is  chosen  from  the  roots  of  a quartic  equation  [6].  In  addition  to 
the  parameters  of  eigenwaves  themselves  in  the  medium  under  study,  there  are  three  parameters 
defining  the  properties  of  the  presented  beams:  9\,  02,  and  kq.  By  setting  these  parameters  in 
various  ways,  one  can  compose  various  localized  fields  with  very  interesting  properties.  Let  us 
illustrate  this  on  the  case  of  localized  fields  in  free  space. 

In  free  space,  the  integral  representations  of  the  localized  fields  under  consideration  and  the 
orthonormal  beams,  presented  in  Ref.  [9],  differ  only  by  the  values  of  the  integration  limits  in 
Eq.  (1).  In  both  cases,  there  are  beams  with  two  independent  polarization  states — Em  and  Ea 
beams. 

Let  us  consider  time-harmonic  fields  Wj  [Eq.  (1)]  with  9i  = 0,  tt/2  < 02  < i r,  and  «o  = 1, 
i.e.,  with  9'  = 9.  These  fields  are  composed  of  plane  waves  propagating  in  a solid  angle  fl  € 
[27T,  47t].  For  the  sake  of  simplicity,  we  assume  that  the  beam  state  function  u = u(9,  p)  reduces 
to  a constant.  A set  of  these  exact  time-harmonic  solutions  of  the  free-space  homogeneous 
Maxwell  equations  consists  of  three  subsets — termed  “storms”,  “whirls”,  and  “tornadoes”  for 
brevity — for  which  time  average  energy  flux  is  identically  zero  at  all  points,  azimuthal,  and 
spiral,  respectively. 

Let  us  first  set  02  = 7T.  Then,  the  fields  under  consideration  are  composed  from  plane 
waves  of  all  possible  propagation  directions,  i.e.,  = Air.  They  are  in  effect  three-dimensional 

standing  waves  with  rather  involved  structures  of  interrelated  electric  and  magnetic  fields.  Beams 
with  s 7^  0 are  essentially  electromagnetic  whirls  with  azimuthal  energy  fluxes.  For  Ea  and  Ba 
electromagnetic  storms,  both  of  which  are  defined  by  the  zonal  spherical  harmonics  (s  = 0), 
the  time  average  Poynting  vector  S is  vanishing  at  all  points.  The  electric  field  vector  E of  Ea 
storms  has  the  only — azimuthal — component,  whereas  the  azimuthal  component  of  the  magnetic 
field  vector  B is  everywhere  zero.  The  opposite  situation  occurs  with  Ba  storms. 

The  spherical  harmonics  with  s ^ 0 define  electromagnetic  whirls  for  which  the  time  average 
Poynting  vector  S has  the  only  nonvanishing — azimuthal — component.  This  component,  as  well 
as  the  energy  densities  we  and  wm  of  the  electric  E and  magnetic  B fields,  is  independent  of 
the  azimuthal  angle  tp.  The  whirls  with  j > s > 1 have  two  major  domains — above  and  below 
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the  equatorial  plane — with  large  energy  fluxes.  The  whirls,  defined  by  the  sectorial  harmonics 
{j  = s > 1),  have  only  one  such  domain,  and  the  energy  flux  peaks  in  the  equatorial  plane. 

Let  us  now  consider  time-harmonic  fields  Wj  [Eq.  (1)]  with  0\  = 0,  7r/2  < 02  < 7T,  and 
Ko  = 1,  i.e.,  with  O'  = 9,  and  2n  < < 4n.  As  before,  we  assume  that  the  beam  state  function 

v = v{9 , p)  reduces  to  a constant.  In  this  case,  the  field  also  is  highly  localized,  but  the  normal 
and  the  radial  components  of  time  average  Poynting’s  vector  S axe  not  vanishing.  As  a result, 
lines  of  energy  flux  become  spiral,  provided  that  s / 0.  We  refer  to  such  specific  localized 
fields  with  spiral  energy  fluxes  as  electromagnetic  tornadoes.  Their  lines  of  energy  flux  closely 
resemble  spirals,  and  as  9i  tends  to  7r,  the  step  of  these  spirals  decreases.  For  the  fields  defined 
by  the  zonal  spherical  harmonics  (s  = 0),  the  lines  of  energy  flux  lie  in  meridional  planes. 

3.  Evolving  Storms,  Whirls,  and  Tornadoes 

Although,  in  many  cases,  the  presented  time-harmonic  solutions  may  provide  satisfactory  models 
of  real  physical  fields,  more  accurate  models  can  be  obtained  by  integrating  these  solutions 
with  respect  to  frequency.  In  particular,  some  solutions  which  describe  quasimonochromatic 
electromagnetic  beams  are  obtained  in  Ref.  [6]. 

Let  us  consider  localized  fields  Wj(r,f)  with  three-dimensional  beam  manifold  = Bx 
[tu_,u;+],  related  with  Wj(r,  t)  [Eq.  (1)]  as 

W|(r,  t)  = ^ J r+  Wj(r,  t)du>,  (3) 

where  Ac o = (cj+  — w_)/2.  In  the  case  of  quasimonochromatic  beams,  Ac o •C  u>.  For  the  beams 
under  consideration,  the  amplitude  function  W (0,p)  is  frequency  independent.  If  the  beam 
state  function  v{0,  p)  also  is  frequency  independent,  or  its  frequency  dependence  is  negligibly 
small,  we  have 

/•2tt  f$2 

Wj(r,<)  = e~tu>t  Jo  dp  Je  p)j0  \p0  (r  • k(0,  p)  - ut)] 

xv(0,p)W(9,p)sm0d0,  (4) 

where  cu  = (w+  + u-)/2  and  po  = Acu/cu.  In  free  space,  this  field  is  composed  of  plane  wave 
packets  radially  moving  with  the  light  velocity  c. 

Therefore,  the  field  under  considerarion  is  essentially  an  evolving  whirl  in  the  neighbourhood 
of  the  point  r = 0.  It  varies  in  intensity  as  different  ” peaks  and  valleys”  reach  the  core  region. 
At  — 7r/Acj  < t < tt/Au,  the  main  peak  passes  through  this  region,  and  the  whirl  reaches 
its  absolute  maximum  intensity  at  t = 0.  At  this  moment,  its  field  structure  is  very  close  to 
the  structure  of  the  corresponding  time-harmonic  whirl.  In  particular,  lines  of  energy  flux  are 
circular  for  both  whirls.  At  — 7t/Acj  < t < 0 and  0 < t < tt/Auj,  the  energy  flux  lines  of  the 
evolving  whirls  are  convergent  and  divergent,  respectively. 

On  the  whole,  the  evolution  of  the  field  can  be  described  as  follows.  When  t — >■  ±oo,  the  field 
tends  to  zero  at  all  points  r.  Therefore,  the  solution  Wj(r,  t)  [Eq.  (4)]  describes  initiation  and 
evolution  of  a whirl,  which  originates  at  the  infinity  at  t = — oo  as  infinitely  small  converging 
wave  propagating  with  the  fight  velocity  c.  At  t -C  — ir/Au,  there  is  a very  small  converging 
wave  with  maximum  peak  at  the  distance  r = —ct.  During  all  this  time,  there  is  also  a weak 
whirl  in  the  neighbourhood  of  the  point  r = 0.  It  passes  through  maximums  and  minimums  of 
activity,  gradually  gaining  in  intensity  as  t tends  to  zero. 

The  total  field  can  be  described  as  the  superposition  of  converging  and  expanding  waves  with 
ever  changing  proportion.  At  t > 0,  the  whirl,  still  passing  through  maximums  and  minimums 
of  activity,  gradually  transforms  into  expanding  wave,  which  vanishes  in  the  infinity  as  t — >■  +oo. 

It  is  significant  that  the  evolving  storms,  whirls,  and  tornadoes  have  finite  total  energy. 
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4.  Conclusion 

In  this  paper,  new  families  of  exact  solutions  of  the  homogeneous  Maxwell  equations  in  complex 
media  and  free  space,  obtained  on  the  basis  of  angular-spectrum  representations,  are  presented. 
These  families  of  solutions  specify  three-dimensionally  localized  electromagnetic  fields  having  a 
rather  involved  curl  structure  and  a very  small — about  several  wave  lengthes  of  composing  plane 
waves — and  clearly  defined  core  region  with  maximum  intensity  of  field  oscillations.  Outside  of 
the  core,  the  intensity  of  oscillations  rapidly  decrease  in  all  directions.  In  complex  media,  these 
fields  provide  a unique  global  description  of  the  medium  under  study,  which  is  supplementary 
to  the  eigenwave  description.  Whereas  each  eigenwave  specifies  the  properties  of  the  medium 
for  one  particular  direction  of  propagation,  the  field  value  of  a three-dimensional  standing  wave 
in  any  point  is  defined  by  all  eigenwaves.  Moreover,  even  in  free  space  such  waves  possess 
very  interesting  properties.  In  a given  Lorentz  frame  L,  a set  of  obtained  time-harmonic  free- 
space  solutions  consists  of  three  subsets  — termed  “storms”,  “whirls”,  and  “tornadoes”  for 
brevity  sake  — for  which  time  average  energy  flux  is  identically  zero  at  all  points,  azimuthal, 
and  spiral,  respectively.  In  any  other  Lorentz  frame  L',  they  will  be  observed  as  a kind  of 
electromagnetic  missiles  moving  without  dispersing  at  speed  V < c.  Solutions  which  describe, 
in  the  frame  L,  finite-energy  quasi-monochromatic  evolving  electromagnetic  storms,  whirls, 
tornadoes,  and  correspondingly,  in  the  frame  L',  various  types  of  moving  and  evolving  missiles, 
are  also  obtained.  The  intrinsic  tensor  technique  in  Minkowski  space  [10]  provides  a convenient 
means  for  investigating  all  these  new  types  of  waves.  Their  properties  are  illustrated  in  graphic 
form. 
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